Drops on an inclined heterogeneous substrate: onset of sliding 
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'U^ ; Abstract 

' Pinning and depinning of drops on an inclined heterogeneous substrate is studied as a function 



of the inchnation and heterogeneity amphtude. Two types of heterogeneity are considered: a 



^ . hydrophobic defect that blocks the droplet in front, and a hydrophilic one that holds it at the 
l> 

lO ' back. Two different types of depinning leading to sliding motion are identified, and the resulting 
O ' 

■ stick-slip motion is studied numerically. 
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It is well known that liquid drops on an ideally smooth substrate move in response 
to external gradients. For example, a drop on an inclined substrate slides downslope in 
response to the gradient of potential energy Likewise a droplet in a temperature 

gradient will move towards higher temperatures as a result of Marangoni forces caused by 
surface tension gradients j^. Alternatively, a wettability gradient induced by a chemical 
grading of the substrate also causes droplet motion. In order to minimize its energy the 
droplet will move towards the most wettable region 0, 0]. Although on ideally smooth 
substrates droplets will move even for arbitrarily small gradients, this is not the case for 
the 'real' substrates used in experiments. There the onset of contact line motion is strongly 
influenced by chemical or physical heterogeneities of the substrate and a finite driving force 
is necessary to overcome the pinning influence of the heterogeneities ^, Q, 0, 0, Ilfl , U , Q ■ 
On the smaller, atomic scale surface heterogeneities can trap droplets even on very smooth 
surfaces. Indeed, heterogeneities occurring on a micro- or mesoscale are known to affect 
the macroscopic movement of droplets and are responsible, for instance, for the observed 
hysteresis between advancing and receding contact angles |l3l. Il4l Il5|. 



The simplest example of depinning is described by the Adler equation 

e = ii-sme, (1) 

where 9 represents the position of the droplet, and fi > represents the applied force. When 
yU < 1 this equation has a pair of fixed points, one of which is stable and the other unstable. 
At /i = 1 these fixed points annihilate in a saddle-node bifurcation, but unlike the standard 
saddle-node bifurcation this bifurcation produces periodic motion for /i > 1. This result is 
simplest to understand if we write Eq. ^ as 6 = —dV/d6, V = —^6 — cos^^. Evidently, 
Eq. represents an overdamped particle in a cosinusoidal potential that is progressively 
tilted as /i increases. A 'particle' in a stable equilibrium at a local minimum of this potential 
'spills out' once the tilt becomes large enough that its position no longer corresponds to a 
minimum. This occurs precisely at /i = 1. The periodic motion present for /i > 1 corresponds 
to the particle sliding down the resulting 'washboard' potential. The period of this motion 
diverges as (/i — 1)"^/^ [17]. The resulting bifurcation is sometimes called a Saddle- Node 
Infinite PERiod bifurcation or 'sniper' for short. 

In this paper we explore the process of pinning and depinning of driven droplets on a 
heteroeeneous substrate. For simDlicitv we consider the case of eravitational forcine on 



an inclined substrate with a heterogeneous disjoining potential with a well-defined spatial 
period such as might arise from spatially varying wetting properties resulting from chem- 
ical heterogeneity. This formulation avoids complications arising from changes in surface 
elevation of the substrate (surface roughness) while retaining the essence of the pinning phe- 
nomenon. In addition we focus on nanoscale droplets for which we can solve the governing 
equation for both the droplet profile and the precursor film, without involving the matched 
asymptotic expansions required for larger droplets. We avoid energy methods since these 
do not permit us to study time-dependent phenomena. 

A two-dimensional liquid droplet on an inhomogeneous solid substrate subject to a hor- 
izontal force /i (Fig^ is described by an evolution equation for the film thickness profile 
h{x,t) derived from the Navier-Stokes equation using the long- wave approximation [l8l |: 

5^/1 = -5, 1^ [5.(75../i + n(/i,x))+/i]|. (2) 

Here 7 is the surface tension, rj is the dynamic viscosity, while Il{h, x) is the disioining 
pressure that accounts for the wetting properties of the heterogeneous substrate We 
use the form Il{h,x) = 2Sadl/ + {Sp{x) / 1) exp[{do — h) / 1] ^Q, 21\, where Sa and Sp are the 
apolar and polar components of the total spreading coefficient S = Sa + Sp, do = 0.158 nm 
is the Born repulsion length and / is a correlation length 2£|, and choose Sa > and 
Sp < 0, thereby combining a stabilizing long-range van der Waals liquid-solid interaction 
with a destabilizing short-range polar interaction. The latter contains the infiuence of surface 
coating and wettability defects, and crucially infiuences the static contact angle When 
fi = the resulting model describes static droplets with a finite mesoscopic equilibrium 
contact angle sitting on an ultrathin precursor film. However, any qualitatively similar 
disjoining pressure yields like results, as shown for dewetting in j2l| and for chemically 
driven running droplets in 22 1. 

We nondimensionalize Eq. using the scales Srj'-f/n'^l for time, I for the film thickness 
and a/TtT^ for the lateral coordinate, where k = ( 15*^1//) exp((io/0- In addition we define 
the dimensionless quantities b = (2Sadl/\Sp\P) exp(—do/l) and a = (7//K^)^/^(/i/p). Thus 
for gravitational forcing a measures the inclination of the substrate, and we refer to it as 
the inclination. The loading of the system (relevant for gravitational forcing) is measured 
by the mean film height h = L^^ h{x) dx. 

Figures Olfa.b) show samDle steadv state profiles for two cases: fa) a hvdroohihc defect. 



(b) a hydrophobic defect. Both are described by 



U{h,x) 



b 



[1 + e^x)] e 



-h 



(3) 



with 



= {2cn[2K{k)x/L,k]Y - A, 



(4) 



where K{k) is the complete eUiptic integral of the first kind and A is such that the average 
of C,{x) over a spatial period vanishes. We use the logarithmic measure s = — log(l — k) 
to quantify the steepness of the heterogeneity profiles (Fig. [2fa,b)). These correspond to 
hydrophilic (e < 0) or hydrophobic (e > 0) defects. In (a) the droplet is held at the back by a 
hydrophilic defect and develops a prominent shoulder as a increases just prior to depinning. 
In contrast in (b) the hydrophobic defect blocks the droplet and its profile steepens with 
increasing a. The profiles at depinning are shaded. Figures |2tc,d) show the advancing and 
receding mesoscopic contact angles, measured at the infiection points of the drop profile, as 
a function of a. For a droplet pinned at the back (Fig. |21^c)) the advancing [receding] angle 
decreases [increases] for small but increasing inclination a. However, once the droplet starts 
developing a shoulder at the back the receding angle decreases again until depinning occurs. 
The situation differs for a droplet pinned at the front (Fig. Efd)). In this case both angles 
increase with a but drop just prior to depinning (Fig.|21). 

The depinning process corresponds to the loss of stability of the pinned drop. The 
stability calculation reveals two mechanisms that lead to depinning. The first is via 
a sniper bifurcation (i.e., a steady state bifurcation) and prevails for hydrophobic defects 
with small wettability contrast 23] and for hydrophilic defects. Figure Efa) shows a typical 
bifurcation diagram for the latter function of increasing a. The figure shows the 

norm of 6h = h{x) — h for pinned drops and its time-average after depinning. Although 
there are two saddle-node bifurcations in the diagram time integration (open circles) shows 
that the upper part of the branch of pinned drops is stable until the rightmost saddle-node 
bifurcation. Thereafter the solutions are time- dependent but periodic (open triangles). The 
inset shows that near the saddle- node the period diverges like (a — a^)"^^^ and hence that 
in this case depinning corresponds to a sniper bifurcation. Figures infb,c) show space-time 
plots of the resulting motion for (b) a > ac, and (c) a = 0.04. In (b) the drop spends a 
long time in a nearly stationary state while slowly spreading downstream, before it abruptly 
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less at constant speed downslope, although the location of the defect remains visible in the 
space-time plot. 

Figure |^a) shows the location of the two saddle-nodes in the (e, a) plane. In the case 
of a hydrophilic defect (e < 0) the saddle-nodes are always present; the one at larger a 
corresponds to the depinning bifurcation. For fixed h and large L the critical a decreases as 
1/L (not shown), as expected on the basis of simple loading ideas. However, the figure also 
shows that something else happens for sufficiently hydrophobic defects. Here the saddle- 
nodes annihilate at e ~ 0.6, and depinning now occurs via a Hopf bifurcation (dashed line). 
The resulting bifurcation diagram (Fig. El^a)) shows that the range of stable pinned profiles 
overlaps with the range of periodic states generated by the instability. Thus in this case 
the branch of periodic solutions loses stability at a saddle-node bifurcation as a decreases, 
and the system settles into a steady pinned state in a hysteretic transition. Figures IStb,c) 
show space-time plots of the periodic state near this transition and further away. Here the 
depinning is as abrupt as in Fig. Efb) but without the slow downstream leakage seen in the 
latter figure. The resulting dynamics strongly resemble stick-slip motion. However, further 
away from the transition the depinned states in both cases look quite alike: in both cases 
the droplet travels at almost constant speed, only slightly modulated by the heterogeneity. 

The advancing and receding angles at depinning (shaded profiles in Fig. Ej) shown in 
Fig.m^b) provide a measure of the contact angle hysteresis observed macroscopically. In the 
case of a hydrophobic defect at the front (e > 0) both angles increase nearly linearly with 
defect strength, and continue to do so even for oscillatory depinning (e > 0.6); the small 
hook visible in the figure near this transition indicates that the Hopf bifurcation sets in 
prior to the disappearance of the saddle-node bifurcations. The behavior is more intricate 
when the pinning is by a hydrophilic defect at the back (e < 0). In this case the role of the 
two angles is reversed, and both decrease nearly linearly with slopes identical to those in 
the e > case. For e < —0.2, however, the receding angle reverses tendency and starts to 
increase again, while the advancing angle continues to decrease. This change in behavior is a 
consequence of the stretching of the drop with increasing inclination just prior to depinning: 
for e < —0.2 gravity drags the main body of liquid downstream (to the right) but the spot 
of higher wettability traps part of it upstream. For fixed a the latter effect becomes more 
pronounced as |e| increases, cf. Fig.^l^b). 

We have examined two types of pinning: pinning by a hydrophilic defect at the back of the 



droplet, and pinning by a hydrophobic defect in front of it, and identified two mechanisms 
whereby pinning takes place. In the case of a sufficiently large hydrophilic defect the droplet 
stretches markedly just prior to depinning as the substrate inclination increases; the inclined 
droplet loses stability at a saddle-node bifurcation, resulting in periodic motion as the droplet 
slides over a periodic array of hydrophilic defects. We have referred to this type of bifurcation 
as the 'sniper'. The periodic motion that results is slow when the droplet is stretching, and 
fast once the droplet breaks away from a defect and spills onto the next one. The situation 
is richer for hydrophobic defects that pin the droplet by blocking it. In this case in addition 
to the steady state sniper bifurcation a new depinning mechanism was observed: the droplet 
loses stability to an oscillatory mode prior to depinning. A mode of this type cannot be 
identified by standard energy arguments. In the example shown this bifurcation is hysteretic. 
The two depinning scenarios are distinguished primarily by the average speed of the droplet 
near the depinning transition. In the sniper scenario this speed vanishes as (a — Oc)^^^; in 
the latter it is finite. At larger values of a both scenarios lead to broadly similar dynamics: 
more-or-less uniform sliding motion modulated by passage over defects. It is noteworthy 
that no Hopf bifurcation occurs when the wettability profile is sinusoidal 2^ . 

Many depinning phenomena in physics may be understood using the sniper scenario. Usu- 
ally this is so in systems with a continuous symmetry such as invariance under translations. 
In the absence of a heterogeneity spatially periodic structures may undergo a spontaneous 
parity-breaking bifurcation that breaks the left-right symmetry of the pattern and produces 
a drift. The direction of the drift is then determined by the associated tilt of the structure 



24| . In this case the drift speed of the structure vanishes as the square root of the distance 



from the parity-breaking bifurcation. However, in the presence of spatial heterogeneities the 
situation changes dramatically because near the bifurcation even small amplitude hetero- 
geneities suffice to pin the tilted structure. A detailed study of this regime shows that 
while some depinning events are indeed analogous to the behavior described by the Adler 
equation, a quite different depinning mechanism is present as well. Here the tilted state 
first undergoes a Hopf bifurcation that produces back-and-forth rocking motion of the tilted 
structure, but no net translation. As a parameter increases the amplitude of this oscillation 
increases leading to a global bifurcation involving an unstable fixed point and its translate 
by one period. This bifurcation generates oscillations with a nonzero mean drift, and this 
net drift increases with further increase in the parameter. The present system differs in the 



absence of left-right symmetry, but a global bifurcation that changes the topology of the 
limit cycle produced in the Hopf bifurcation from a libration to a rotation must still take 
place. Such a bifurcation can occur if the Hopf bifurcation is in fact supercritical. Consis- 
tency with the Fig. Efa) requires that the branch of periodic states must go through a pair 
of saddle-node bifurcations to produce stable states of the type shown in Fig. E^b). 
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FIG. 1: Sketch of a droplet on a heterogeneous substrate subject to a horizontal force /i towards 
the right. 




FIG. 2: Characteristics of pinned droplets as a function of the forcing a for locahzed hydrophihc 
[(a) and (c), e = — 1] and hydrophobic [(b) and (d), e = 1] defects. The upper parts of (a) and (b) 
show steady droplet profiles while the lower parts show the wettability profile [Eq. (jlj with s = 6]. 
The profile at depinning is shaded. In (a) the droplet is pinned by a more wettable defect at the 
back whereas in (b) it is blocked by a less wettable defect in front. Panels (c) and (d) show the 
advancing and receding contact angles ^ as a function of a. The remaining parameters are L = 25, 
b = 0.1, and h = 1.5. 
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FIG. 3: (a) Bifurcation diagram for depinning via a sniper bifurcation for a hydrophilic defect 
[Eq. @ with s = 6] with e = —1.0 and L = 25, h = 0.1, h = 1.5. The figure shows the L^-norm of 
steady solutions (sohd hne), selected steady solutions as obtained by integration in time (circles) 
and the time-averaged L2-norm for the unsteady solutions beyond depinning (triangles). Inset 
shows the inverse of the temporal period T for the latter. The remaining panels show space-time 
plots over one spatial and temporal period for a sliding drop (b) close to depinning at a = 0.0185 
with T = 556.1, and (c) far from depinning at q = 0.04 with T = 100.7. 
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FIG. 4: (a) Phase diagram for the depinning transition for localized hydrophilic (e < 0) or hy- 
drophobic (e > 0) defects [Eq. ^ with s = 6]. The figure focuses on small wettability contrast and 
L = 25, b = 0.1, h = 1.5. The solid [dashed] lines correspond to saddle-node [Hopf] bifurcations. 
The latter emerge near the cusp at which the two saddle-node bifurcations annihilate for e > 0. 
(b) Advancing (solid lines) and receding (broken lines) contact angles 9 at the depinning transition 
as a function of wettability contrast for a hydrophilic defect at the back (e < 0) and hydrophobic 
defect at the front (e > 0). Thick [thin] lines refer to depinning through a real [oscillatory] mode. 




FIG. 5: As for Fig. |21but showing depinning via a Hopf bifurcation when e = 1.0. (a) Bifurcation 
diagram, (b) Space-time plot for a = 0.0415 with T = 206.4. (c) a = 0.08 with T = 47.4. The 
vertical line indicates the location of the Hopf bifurcation as obtained from Hnear stability theory 



